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Abstract 

Graphene nanoribbons (GNR) in mutually perpendicular electric and 
magnetic fields are shown to exhibit dramatic changes in their band struc- 
ture and electron transport properties. A strong electric field across the 
ribbon induces multiple chiral Dirac points, closing the semiconducting 
gap in armchair GNR's. A perpendicular magnetic field induces partially 
formed Landau levels as well as dispersive surface-bound states. Each 
of the applied fields on its own preserves the even symmetry Ek = E-k 
of the subband dispersion. When applied together, they reverse the dis- 
persion parity to be odd and gives Ee,k — —Eh,-k and mix the electron 
and hole subbands within the energy range corresponding to the change 
in potential across the ribbon. This leads to oscillations of the ballistic 
conductance within this energy range. 

Recent advantages in the fabrication techniques of graphene nanoribbons 
(GNR) together with the long electron mean free path have stimulated consid- 
erable interest in their potential applications as interconnects in nano circuits. 
Near the K and Dirac points for infinite graphene, the electrons are massless 



and chiral |Neto et al.(2009)Neto, Guinea, Peres, Novoselov, and Geim . The 



electronic properties of GNR are sensitive to the geometry of their edges and the 
number of carbon atoms N across the ribbon. The GNR is thus classified as arm- 
chair (ANR), zigzag (ZNR) nanoribbons for even N and their counterpart anti- 
armchair (AANR), anti-zigzag (AZNR) for odd TV. The armchair confinement 
mixes K and K' valleys creating chiral electrons around the F point. Chirality 



is the key ingredient for unimpeded electron transport (Klein effect). Depend- 
ing on modulo 3, the ANR/AANR can be either metallic or semiconducting 
making them suitable candidates for use as field- effect transistors. In contrast, 
the zigzag confinement does not mix the valleys but rather intertwine their longi- 
tudinal and transverse momenta, creating edges-bound quasi-particles between 
the K and points. For ZNR/AZNR, the electrons are not chiral (in the sense 
of projection of the pseudo-parity on the particle momentum), and the electron 
transmission through a potential barrier is determined by the electron pseudo- 



parity Rainis et al.(2QQ9)Rainis, Taddei, Dolcini, Polini, and Fazio] . This quan- 



tity redefines the Klein effect as the suppressed transmission through the barrier 



in ZNR, also known as the valley- valve effect Roslyak et al.(2009)Roslyak, lurov, Gumbs, and Huang 
The latter is the basis for the proposed valley filters. The electron confine- 
ment in GNR causes their properties to be quite sensitive to an applied electric 



Novikov(2QQ7)| [Novikov and Levitov(2QQ6)| |Raza and Kan(2QQ8)| or magnetic 



Brey and Fertig(2QQ6)[|Ritter et al.(2QQ8)Ritter, Makler, and Latge[|Perfetto et al.(2QQ7)Perfetto, Gonzalez, ( 



Golizadeh-Mojarad et al.(2QQ8)Golizadeh-Mojarad, Zainuddin, Klimeck, and Datta] 



field. These changes are reflected in measurable quantities such as the bal- 
listic conductivity and local density of states (LDOS) [Lyo and Huang (2004) [ 



Li et al.(2009)Li, Huang, Chang, Chang, and Lin|. 



In this letter, we report on the individual and combined effects of an electric 
and magnetic field on the band structure and conductance of GNRs. If only one 
of the fields is applied, it is well known that the time reversal symmetry of the 
energy bands for electrons and holes is preserved for all the types of GNRs we 
listed above. However, the combined effect of an electric and magnetic field on 
the energy is to break the time reversal symmetry for both electrons and holes 
and mix the energy bands. The effect of mixing on the differential conductance 
and LDOS is presented below and our results are compared to those obtained 
when only one of the two external fields is applied to an ANR with quantum 
point contacts as illustrated schematically in Fig. [l] The ribbon is attached to 
left (L) and right (R) leads serving as infinite electron reservoirs. The R-lead 
is assumed to be the drain held at chemical potential ja. The L-lead is held 
at DC biased chemical potential /i + eV (e is the electron charge and V is the 
bias potential) and serves as the source. We choose coordinate axes so that the 
nanoribbon is along the x axis in the x?/-plane. Mutually perpendicular static 
electric field Ey along the y axis and magnetic field Bz along the z axis are 
applied, as shown in Fig. [ij 

We calculated the energy bands for graphene with sublattices A and B in the 
tight-binding model [Lin et al.(2009)Lin, Hikihara, Jeng, Huang, Mou, and Hu] 



Neto et al.(2009)Neto, Guinea, Peres, Novoselov, and Geim|. These are then 



separated into hole {h} = {1 < n < N} and electron {e} = {A^ < n < 2N} 
energy bands. The two component wave function is a normalized 2N vector 
{'^{x)\n^k = ((^A(^)|n,/c7 (^B(^)|n,/c)- The clcctric field induces a potential 
across the ribbon U{y) = —eSy{y — W/2) = Uo{y — W/2)^ where W is the ribbon 

■"^ Since we neglect spin, the action of the time reversal operator T amounts to reversing 
the direction of the wave vector propagation. The even/odd particle energy symmetry may 
be defined as En^k = ^TEn^k = ^En^-k 
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Figure 1: (Color online) Schematic of an ANR in the presence of an in-plane 
electric field E along the y axis and a perpendicular magnetic field B along the 
z axis. 

width. The magnetic field modifies the wave vector diS k ^ k — BzCy /c, which 
amounts to the Peierls phase in the hopping integrals |Liu et al. (2008) Liu, Ma, Wright, and Zhang . 
The magnetic field strength is assumed weak so we could take the energy lev- 
els as spin degenerate. The dispersion curves can be experimentally observed 
via scanning tunneling microscopy (citation). The tunneling current flowing 
through the microscope tip is proportional to the LDOS given by 




(1) 



n,k 



The energy dispersion determines the ballistic charge transport through the 
ribbon, at temperature T, by 



2e 



Q{V,fi,£y,B,,T) = 




(2) 



where Vn,k = dEn^k/d{hk) is the carrier group velocity. At T = 0, the Fermi 
function at the source contact is /^^ = 1 — 6> {En^k — M — eV) and for the drain, 
we have /^^ = 1 — {E^^k — 1^)- We note that Eq. ([2| does not assume any 
symmetry for the energy dispersion relation. If the energy satisfies E^^k = 
En^-k^ we obtain the well-known Landauer-Butikker formula [Buttiker(1988)| . 
That is the differential conductance G{ii^£y^Bz) = {dl /dV)y^Q is determined 



3 



by the number of right-moving carriers through i'n,/c/|'^n,/c| > at the chemical 
potential E^^k = Alternatively, one may take the difference between the local 
minima and maxima below the chemical potential En^k < M (citations). 

Our numerical results for the energy bands, LDOS and conductance for semi- 
conducting ANR {N = 51) in the presence of an electric and/or magnetic field 
are presented in Figj2j When either only an electric or magnetic field is applied 
£yBz = 0, the electron/hole energy bands are symmetric with E^^j. = — ^e,/c 
and time reversal symmetry is satisfied with En^k = En-k around the /c = 
Dirac point in Figj2|b.l). The latter means that if the time for the particle is 
reversed, the particle retraces its path along the same electron/hole branch. The 
LDOS also demonstrates the wave function symmetry with respect to the rib- 
bon center LDOS(£:,x) = LDOS(-£;,x) = LDOS(£:,-x). In accordance with 
the Landauer-Butikker formalism, the conductivity demonstrates the familiar 
staggering pattern. The magnetic field by itself distorts the weak dispersion 
(n close to N) so that the partially formed Landau levels E^^q ~ \/B/ri shows 
itself up as the fiat parts in the dispersion curves. The lowest Landau level 
provides the single conducting channel (along the ribbon edges), while the rest 
are doubly degenerate. When the wave vector evolves from the Dirac point, 
the degeneracy is lifted and the lowest subband acquires a local minimum. Of 
these two effects, the first one can be observed in the LDOS, while the second 
reveals itself as sharp spikes in the conductance as depicted in Fig. [2|b.3). For 
the high energy subbands, when the radii of the Landau orbits (spread of the 
wave function in Fig. |2] (b.2)) become comparable with the ribbon width, the 
confinement effects dominate and the spectra become linear in magnetic field 
with En^Q ^ Bz/n. These subbands are not degenerate. 

The main effect which the electric field has on the energy dispersion is to frac- 
ture Fermi surface into small pockets for /c 7^ 0, and thereby closing the semicon- 
ducting energy gap. These zero energy points, where the group velocity abruptly 
changes sign, represent new Dirac points, which follows from the chirality of the 
wave function in their vicinity Brey and Fertig(2009)t . The rapid changes in the 
group velocity cause the appearance of spikes in the conductance near < Uq/2 
and its step-like pattern is broken. Due to the Dirac symmetry of the problem, 
the electron-hole band structure remains symmetric. The energy dispersion is 
not affected by magnetic field at the original Dirac point /c = 0. Time reversal 
symmetry also persists. The LDOS shows that at high energies the electric field 
confines the electrons and holes near opposite boundaries. However, at low ener- 
gies the LDOS does not change across the ribbon, which is a manifestation of the 
Zitterbewegung effect (attempt to confine Dirac fermions causes wave function 
delocalization Neto et al.(2QQ9)Neto, Guinea, Peres, Novoselov, and Geim| ). With 
respect to the three cases considered above, we point out that the hallmark of 
Dirac fermions is the even symmetry of the dispersion with respect to the wave 
vector, and steams from time reversal symmetry. Even though an attempt to 



confine them may lead to the broken electron/hole symmetry Peres et al.(2006)Peres, Castro Neto, and Guinea 
the wave vector symmetry still persists. 

We now turn our attention to the most interesting case when both elec- 
tric and magnetic fields are applied together. Concurrent action of the electric 
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field dragging force, tiie Lorentz force and confinement by the ribbon edges de- 
stroys the Dirac symmetry of tiie problem so that £^n,/c 7^ ^n,-/c ^is shown in 
Figj2|d.l). The dispersion distortion is different for the electrons and holes, so 
the symmetry between the conduction and valence bands is also broken. On 
one hand, the partially formed Landau levels get distorted by the confinement 
due to the electric field in conjunction with the edges. Their degeneracy is also 
lifted. On the other hand, the magnetic field does not allow formation of addi- 
tional Dirac points and wave function delocalization. At high energies, where 
the group velocity is decreased and the drag due to the electric field prevails. 
The electrons and holes get gathered at the opposite ribbon edges (Figj2|d.2)). 
For lower energies, in the region \En^k\ ^ ^o/2, the electron/hole dispersions 
overlap. The electron bands have only local minima, whereas only the hole 
bands have local maxima. Regardless of the broken Dirac k symmetry of the 
dispersion, our numerical simulation of the differential conductivity shows that 
the Landauer-Butikker expression still applies. Therefore, in the overlapping 
region < /7o/2, the conductivity oscillates since the minimum of the electron 
band is followed by the maximum on the hole band when the chemical potential 
grows. As for possible applications of the broken Dirac symmetry, the ribbon 
, subjected to mutually transverse electric and magnetic fields, may serve as a 
field-effect transistor with a tunable working point. An interesting feature of 
our results is that there is not only a breakdown in the even-/c symmetry of the 
energy dispersion relation, but the energy bands are reversed with odd symme- 
try, satisfying E^^k = —Eh-k- We explain this effect by adopting the method 
described in Ref. Novikov and Levitov(2QQ6)[ |Brey and Fertig(2QQ9)| . Let us 
focus on the energy region close to the original Dirac point = 0, where the 
unperturbed wave functions are governed by the conventional Dirac equation. 
Both applied fields are treated perturbatively. The effect of magnetic field is in- 
cluded through the wave vector replacement k ^ k — BzCy/c. The electric field 
is treated by a chiral gauge transformation. This transformation shows that the 
spectrum at /c = is affected by the electric field in the presence of the magnetic 
field. Regardless of the metallic or semiconducting ANR the electron and hole 
dispersion become degenerate around k = with Ee^k = —Eh-k ^ — {^y/^z) k. 

In conclusion, we have demonstrated that when GNRs are placed in mu- 
tually perpendicular electric and magnetic fields, there are dramatic changes 
in their band structure and transport properties. The electric field across the 
ribbon induces multiple chiral Dirac points, whereas a perpendicular magnetic 
field induces partially formed Landau levels accompanied by dispersive surface- 
bound states. Each of the fields by itself preserves the original even parity of 
the subband dispersion, i.e. En^k = E^-k^ maintaining the Dirac fermion sym- 
metry. When applied together, their combined effect is to reverse the dispersion 
parity to being odd with Ee^k = —Eh-k and to mix electron and hole subbands 
within an energy range equal to the potential drop across the ribbon. Broken 
Dirac symmetry suppresses the wave function delocalization and the Zitterbe- 
wegung effect. The Butikker formula for the conductance holds true for the odd 
k symmetry. This, in turn, causes the ballistic conductance to oscillate within 
this region which can be used to design tunable field-effect transistors. 
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k #of carbon aloms across AN R dl/dV(V=0) 



Figure 2: (Color online). Panels (a) represent the dispersion curves for the 
electrons (green curves) in the conduction band and holes (blue curves) in the 
valence band. The lowest conduction and highest valence subbands are given 
by the thick curves. Panels (b) show local density of states. Panels (c) give the 
corresponding ballistic conductance in units of 2e^ /h. Panels (1) correspond to 
absence of the em. field. Panels (2) correspond to the sole magnetic field of the 
fiux through a single hexagon placket (/)/0o = 1/150. Panel (3) show the effect 
of the sole electric field of the strength Uq/Iq = 1/2. Panels (4) demonstrate the 
combined effects due to an electric and magnetic field with the same strength 
as that employed in panels (2)-(3). g 
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